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ABSTRACT 

Context. The characterization of a time series is a powerful tool for investigating the nature of mechanisms that generate variability in 
astrophysical objects. Blazar variability across the entire electromagnetic spectrum is a long-standing puzzle, and it has been difficult 
to ascertain the mechanisms at play. 

Aims. Lognormal variability in X-ray light curves, probably related to accretion disk activity, has been discovered in various compact 
systems, such as Seyfert galaxies and X-ray binaries. Identifying a similar behaviour in blazars would establish a link between them. 
Methods. Public X-ray data from the blazar BL Lac are used to investigate the nature of its variability, and more precisely the flux 
dependency of the variability and the distribution of fluxes. 

Results. The variations in the flux are found to have a lognormal distribution and the average amplitude of variability is proportional 
to the flux level. 

Conclusions. BL Lac is the first blazar in which lognormal X-ray variability is clearly detected. The light curve is orders of magnitude 
less variable than other blazars, with few bursting episodes. If this defines a specific state of the source, then the lognormality might 
be the imprint of the accretion disk on the jet, linking for the first time accretion and jet properties in a blazar. 
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1. Introduction 

The luminosity fluctuations from astronomical sources are of- 
ten used as a probe to constrain the physical mechanisms that 
generate them. Nonthermal X-rays from blazars are understood 
to originate in relativistic electrons cooling either through the 
synchrotron or the inverse Compton mechanisms. Their spec- 
tra are usually featureless power law continuums, with little 
or no thermal emission from the disk surrounding the cen- 
tral black hole. If thermal emission is present, it is difficult 
to find because of the dominance of the nonthermal radiation 
( |Perlman, Addison et al. 2008] >, complicating further the study 
of the link between accretion and jet physics in these objects. 
Analyses have attempted to identify in particular features in light 
curves where the flux varies significantly, such as peaks, inter- 
preted as the light crossing time through an emission zone and 
hence a direct probe of its geometrical size. The cause of vari- 
ability is still poorly understood despite decades of observations 
across the entire electromagnetic spectrum. 

We analyze the light curve of BL Lacertae (BL Lac), a mod- 
erately X-ray bright blazar of the class named after itself. It 
has never shown dramatic X-ray outbursts, where fluxes vary 
by over an order of magnitude or more, as is often the case for 
this class of source when they are also very high energy y-ray 
emitter^]. The available data set has already provided a wealth 
of results, the most notable being the observation of a double X- 
ray flare in 2005, interpreted by Mars cher et al. (2008)| as origi- 
nating from a disturbance passing through two different zones in 
the jet of BL Lac. The relatively major outburst-free light curve 
is probed for evidence of a statistical property found in some 
accreting sources, called lognormality, found in galactic as well 
as extragalactic accreting sources which can be generated by a 
stationary process by taking the exponential of a Gaussian time 



BL Lac's spectrum extends up to ~ 1 TeV, see Albert et al. 2007 



series ( |Uttley et al. 2005[ ). However, [McHard y (2008) suggested 
that variations generated within the accretion disk can modulate 
the nonthermal jet emission. Evidence of lognormal variability 
has also been found in very high-energy gamma-ray emission 
from the BL Lac object PKS 2155-304 ( [Degrange et al. 2008) , 
Lognormal fluxes have fluctuations that are on average propor- 
tional to the flux itself, and are the signature of an underlying 
multiplicative physical process, rather than additive. This signa- 
ture is particularly difficult to find in blazars because they are 
relatively faint, causing the Poisson noise to be large, and usu- 
ally densely sampled only on flaring occasions and for a few 
days, but even then fluxes remain modest compared to those of 
X-ray binaries. 

2. Data and analysis 

BL Lac has been observed routinely by the PCA detectors 
dJahoda et al. 1 995) onboard RXTE for many years, and we an- 
alyze 6.75 x 10 5 s of public data collected since March 2005, 
with exposures of a few ks per observation. This corresponds 
to the RXTE ObsID's numbered 91127, 92103 and 93135. The 
PCA data were analyzed using the HEASOFT 6.5.1 package ac- 
cording to the guest observer facility recommended criteria, and 
limited to the unit PCU2 because of known background control 
issues with unit PCU0. Although the RXTE monitoring of BL 
Lac did not occur on a daily basis, a light curve is derived in 1 
day bins (top panel in Fig.Q}, yielding 402 flux estimates fa (and 
associated measurement errors cr,). In order to apply Gaussian 
statistics to our variance measurements, this light curve is di- 
vided further into segments of equal length such that each con- 
tain at least N = 20 points (middle panel) of arithmetic mean 
O. The shortest timescale allowing this condition is found to be 
T — 77 d, where most of the segments have N > 20 flux points. 
The excess rms cr X s, defined as o~ xs = jj Yj'^iKfa ~ ®) 2 ~ of] 
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Fig. 1. Top panel: PCA flux light curve (in units of counts s _1 ) of 
BL Lac with bins of 1 day. Middle and bottom panel: PCA mean 
flux and excess rms measured from 77 day segments (see text). 
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Fig. 2. Distribution of the fluxes for the considered period (top) 
and the distribution of the logarithm of the fluxes (bottom). The 
lines are the result of a Gaussian fit to the data. The improve- 
ment to the fit provided by the logarithmic function is clearly 
apparent in that it reproduces more successfully the steeper-than- 
Gaussian left-hand part of the distribution, as well as the broader, 
tail-like right hand part, which are characteristic of the lognor- 
mal distribution. The single flux estimation of the flare on MJD 
54075 has been excluded so as not to favour the logarithmic fit. 



( |Vaughan et al. 2003] ), which estimates the rms corrected for the 
Poisson noise, is then derived for each of the segments (bottom 
panel). Varying this length T, or requiring a larger amount of 
points N in the variance estimation, does not change significantly 
the results presented here. The error in the excess variance is the 
68% confidence level estimated from 10 4 simulated light curves, 
where a random flux from a Gaussian of dispersion equal to the 
flux error <x, is added to each individual flux estimate <p t (as in 



Appendix B in Vaughan et al. 2003 1. Its relatively flat light curve 
compared to other X-ray blazars, along with the highest flux 
point visible at MJD 53075, where the flux more than tripled, 
indicates that the source spends most of its time in a nonflaring 
state, which we refer to as its "quiescent" state. The short fore- 
mentioned flare is clearly an exception in the overall variability 
of this source. A more detailed analysis of this single 2.4 ks ob- 
servation shows that the 2-10 keV flux of 5.2xl0~ n ergcirT 2 s _I 
was associated with a spectral power law index of F = 2.5 + 0.1, 
much brighter than the 1 -3x 10~ n ergcm 2 s _I fluxes reported 
in Fig. 2 of Marscher et al. (2008), but also very soft compared 
to its brightness. Since our aim is to characterize the long-term 
lightcurve, this single observation is excluded from the analysis, 
removing one flux estimation out of 402. It becomes obvious 



why removing this event is also conservative in terms of estab- 
lishing the nature of the variability. 

The light curve fluxes are presented in a histogram in Fig. [2] 
(top panel), along with the result of a chi-square Gaussian fit. 
The fit of the distribution of the fluxes has a chi-square proba- 
bility of p(x 2 ) = 0.015 (x 2 = 26.4, n.d.f = 13), and a width of 
cr = 0.26 s , indicating that statistical fluctuations alone can- 
not explain the fit. The natural logarithm of the fluxes provide 
a much better Gaussian fit with p(x 2 ) = 0.47 (y 2 = 12.75, 
n.d.f. = 13) (bottom panel), and a width of cr — 0.26. The scatter 
plot of <Jxs versus the average flux in Fig.[3j in which the short 
flare is excluded, shows a clear correlation (p = 0.73). A fit to a 
lineOt 2 = 33, n.d.f = 12) yields crj^ cc (0.15+0.02)0. An F-test 
shows that the improvement to a line fit is significant at a confi- 
dence level exceeding 99% compared to a constant fit (x 2 = 85, 
n.d.f = 13). Including the flare observation causes the linear fit 
to be further preferred to the constant fit, and increases the slope 
of the linear fit by ~ 4%. It is important to remark that the flare 
observation is statistically very unlikely to be just a large varia- 
tion of the same lognormal process, which is another indication 
that it is more likely to be a physically distinct state. 
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3. Discussion 

The lognormality of the flux distributions, and the evidence that 
the fluctuations in the flux are proportional to the flux, imply 
that the variations are lognormal ( Aitchin son et al. 19 63). The 
proportionality factor of 0.15 in Fig. [3] is comparable to that 
(0.17 with a 1-day mean flux binning using ROSAT data, and 
0.26 with a 0.5-day binning using ASCA data) found in the most 
X- ray variable nar row-line Seyfert 1 galaxy IRAS 13244-3809 
by Gaskell (2004), the only other known case of an extragalactic 
source exhibiting lognormal X-ray flux variations. The width in 
Fig. [2] of cr = 0.26 is however smaller (0.44 dex and 0.54 dex 
for the two different data sets mentioned above) than in IRAS 
13244-3809, meaning that the fluctuations are smaller in this 
BL Lac than in the Seyfert galaxy. A linear correlation be- 
tween the excess rms and the average flux was also found in 
the Seyfert 1 galaxy Mrk 766 ( Vaughan et al. 2003 1. The physi- 
cal process responsible for X-ray emission in Seyfert 1 galaxies 
is probably thermal emission from the accretion disk. We also 
note that a linear flux-rms relationship was found for Galactic 
accreting sources such as X-ray binaries Cyg X-l and SAX 
J 1808.4-3658 by |Uttley & M c Hardy (2001)) of proportional- 
ity factors of 0.33 and 0.31, respectively. A specific observa- 
tion of Cyg X-l exhibited a remarkable lognormal fit to its X- 



ray flux distribution ( Uttley et al. 2005)>. The origin of th e linear 
flux-rms was interpreted by |Uttley & M c Hardy (2001) as orig- 
inating from the subdivision of magnetic reconnection energy 
release as an avalanche occuring on large scales in the corona. 
The same authors note that the radius-dependent fluctuations in 
the mass accretion rate in the model of Lyubarskii et al. (1997)] 
can also explain this relationship. It seems possible that the 
photon breeding mechanism dStern & Poutanen 2008 , and ref- 
erences therein), which generates highly non-linear light curves 
in relativistic jets and operates in sites that are not necessarily 
directly connected to accretion properties, could be an alterna- 
tive explanation of lognormal blazar variability, but this remains 
to be investigated since the statistical properties of the photon 
luminosity from this model are not known at this point. 

These striking similarities suggest that the variability mech- 
anism causing the X-ray modulation in BL Lac, outside of 
flaring events, is independent of the emission mechanism, 
and originates in the disk, rather than in the jet, provided 
these lognormal fluctuations are indeed present in the disk. 
The nonthermal emission from BL Lac is usually modeled 
by the synchrotron self-Compton scenario (SSC, see e.g., 
IMaraschi, Ghisellini & Celotti 1992| |Bicknell & Wagner 2 002 ), 
where relativistic electrons cool either via the synchrotron or 
the inverse Compton mechanism. The comptonized photons are 
assumed to be predominantly the synchrotron photons them- 
selves, rather than photons coming from the outside the emis- 
sion region as in the external Compton model (EC, see e.g. 
Sikora et al. 1994| l, given the absence of thermal signatures in 
the emission spectrum of BL Lac. However, these mechanisms 
do not provide a way to generate lognormal variability if one of 
the parameters affecting the predicted flux (e.g., magnetic field, 
Doppler factor, injection rate) does not have a lognormal time- 
dependence itself. If fluctuations in the disk are indeed the origin 
of the lognormal radiation in BL Lac, then the most appropriate 
hypothesis would be one in which lognormal fluctuations in the 
accretion rate yield an injection rate with similar properties. If 
however external photons are at least partly responsible for the 
cooling of the X-ray emitting electrons, then lognormal fluctu- 
ations in the disk luminosity would have a similar result. An 
ambitious observation programme targeting other X-ray bright 
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Fig. 3. Scatter plot of the excess variance versus the average of 
the fluxes for which the excess was determined. The line is a 
linear fit showing cr xs oc (0.15 + 0.02) <I>. 



blazars, over similar time scales, and with a similar sampling 
quality, would allow us to search for lognormal variability, and 
perhaps solve an important piece in the puzzle about blazar vari- 
ability. 
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